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Abstract 

In the framework of the approximation of slowly varying amplitudes 
a multiwave dynamical theory of neutron diffraction on a moving phase 
grating was developed. The influence of the velocity of the grating, its 
period and height of the slits on the discrete energy spectrum and intensity 
of various diffraction orders was analyzed. 


1 Introduction 


In [l] it was predicted that a diffraction grating moving across the beam 
of slow neutrons, is a time-dependent quantum modulator that converts the 
energy spectrum of incident neutrons in a series of quantum states with discrete 
energies. Later a discrete spectrum in the diffraction of ultracold neutrons 
(UCN) on a moving grating was observed in the experiments [2j|5]. In addition, 
it was shown that a moving grating with a varying spatial period can serve as 
a quantum time lens for UCN |6 . 

Let us briefly present the result of the corresponding quantum problem. In 
the laboratory coordinate system a plane neutron wave H h n (x, z , t ) = Aoexp(ikoZ- 
iwot) is assumed to propagate along the z-axis towards a periodic grating, which 
is supposed to be normal to the neutron propagation axis (Fig. |la| ), where 
ko = mVo/H is the neutron wave number, m is the neutron mass, h is the 
Planck constant, loq = hk^/Zm. The grating moves with a constant velocity V 
in the positive direction along the x axis (Fig. lb) and the grating grooves are 
oriented along the y axis. As in 


[lj|4j let us solve the problem in the moving 
coordinate system (a/, z) where the grating is at rest. In this system the wave 
is obliquely incident on the grating: 


% n {x', z , t) = A in exp(ik 0 z - iw't), 


(1) 


with Ai n (x) = Aoexp(ikvx), ky = mV/h , uj' = to o + wy, uy = hky/2m. After 
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the passage of the grating in the region z >0 we have 


OO 

^'(x', z, t) = Aq J F(q)exp(iqx' + ik z z — iui't) dq, 


( 2 ) 


where k z = (/cq + ky — q 2 ) 1 ^ 2 , F(q) is the Fourier transform of the product 
F(x) = A in { x)f{x), where f(x) is the grating transmission function. For a 
periodical function f{x) = JA djexp(iqjx) where j are integers, qj = 2tt j/d, 
d is the spatial period of the grating, F(q) = A 0 JA <5(<? — qj + ky), and the 
Fourier coefficients are 


d 

aj = 2 J f{x)exp(iqjx) dx, j = 0, ±1, ±2 ,... (3) 

o 

Returning to the laboratory coordinate system x = x' + Vt we obtain a super¬ 
position of plane waves with amplitudes cij , discrete frequencies u>j and wave 
vectors ky = ( qj,k Z j) (see Fig, [lc]): 

y(x,z,t) = E a,jexp(iqjX + ik z jZ — iujt), (4) 

3 

where k Z j = (fc§ + 2 kyqj — qj) 1 ^ 2 , u)j = too + = 2tt/T, T = d/V. We 

should emphasize that the waves with different j do not affect each other. 




Figure 1: Schemes of diffraction, a grating at rest ,tOj = wo; b coordinate system 
moves together with the grating tOj = u/; c moving grating in the laboratory 
coordinate system toj = loq + jfh 


Let us consider a grating with a rectangular groove profile with width d/2 
and depth h, so that the phase difference </> = /co(l — n)h = 7 r, where n is the 
neutron refraction index. At normal incidence the amplitudes of even orders 
(including the zero order) are ay = 0, but odd-order ones are cij = 2i/irj (j = 
2s — 1). It is obvious that for larger values of spectral splitting = 2irV/d 
it is necessary to increase the grating velocity V and (or) decrease the spatial 
period d. But due to the oblique incidence in the moving coordinate system the 
phase profile varies with increasing velocity and takes the form of a trapezium 
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I] (Fig# It follows from © that amplitudes of even orders are aj = cBj and 
those of odd orders are aj = — Bj/j , where 

Bj = [1 + exp{-injc)\/[in(l - j 2 c 2 )], (5) 

where dimensionless parameter c = 2hV/dVo is introduced [3j. At c = 1 the 
phase profile takes the form of a triangle. In the case of a moving grating (c > 0) 
the amplitudes of zero and all even orders never vanish and their intensities 
\aj\ 2 increase with increasing c, contrary to the intensities of odd orders, which 
decrease (see Fig. 3). In the experiment [5] a decrease in the intensity of minus 
fsi orc [ er w jth increasing velocity of the grating was really observed, but it was 
greater than that predicted by the relation ©• 

Quantum modulation of neutron waves by a moving grating is a method¬ 
ological basis of experiments to test the equivalence principle for the neutron 
[5,7. Therefore a more rigorous theoretical study of the phenomenon is re¬ 
quired. In this paper we developed a dynamical theory of neutron diffraction 
on a moving grating which to some extent satisfies this requirement. It is based 
on the method of slowly varying amplitudes and takes into account the mu¬ 
tual influence of waves of different orders as the neutron penetrates a 3D phase 
grating. 


2 Theory 

Neutron propagation in a medium is described by the Schrdinger equation 

A>F(r,t) + [k 2 - fc^(r)]\l/(r,t) = 0, (6) 

where A is the Laplace operator, k is the wave number of neutrons in vacuum, 
k 2 ( r) = 47rA(r)6(r), N{ r) is the density of nuclei, 6(r) is the coherent scattering 
length of the medium. Let us assume that function k 2 (x) periodic in region 
0 < z < h has the form 


where q n 


k%(x) = ^ Xnexp{iq n x), 


27 m/d are reciprocal lattice vectors, 


Xn = 


k 2 {x)exp{—iq n x) d.T. 


(7) 


( 8 ) 
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A zero Fourier amplitude %o determines the value of the average grating 
refractive index n e = (1 — Xo/k 2 ) 1 ^ 2 in a layer of thickness h. 

The wave function of neutrons in the region 0 < z < h in the moving 
coordinate system we can write using ([I]), as the sum of the Bloch functions 
with amplitudes 4/ m (z) depending on the vertical coordinate z: 

OO 

^'{x',z,t)= ^2 ^ m (z)exp{iq m r' - iu't), (9) 


where the projections of the wave vectors q mx = Qm — ky, q m z = k$ z = (/c(j — 
X'o) 1 ^ 2 - Here, it was assumed that the wave number in vacuum is k' = (k$ + 
ky) 1 / 2 , and in the layer 0 < z < h it is k 0z = ( k l2 n 2 — ky) 1 / 2 . 

By inserting Q and into ([6]), we equate the terms with equal exponents 
and neglect the second derivative of ’F m with respect to z. As a result, we 
obtain a system of differential equations: 


d’F„ 

dz 


— i ^ ( f3 n ty m — r 

o 


( 10 ) 


where a m = q rn (2ky — q rn )/2ko Zl /3 n = Xn/2kQ Z . System of equations ( |10| should 
be supplemented by boundary conditions : \Ho(z = 0) = AO, = 0) = 0. 

Neglecting second derivatives when deriving ( |To) ), we excluded the reflected 
waves from consideration. Apparently, this should not result in serious errors, 
since for typical experimental conditions [2j|6| the Fresnel reflection coefficient 
is rather small: \Rf\ 2 « 2 — 4%. 


It can seen from (10) that unlike the simple model described above, the 


amplitudes of waves of arbitrary m-th order depend on the amplitudes of other 
orders. Coupling coefficients f3 n in addition to the order number are also defined 
by material and shape of the grating profile. In our case Xn = k 2 sin(nn/2)/Trn. 
The a m values depend on the grating velocity and the diffraction order. They 
vanish both at m = 0 and Bragg condition 2 ky = q m - The system of equa¬ 


tions (101 is solved using the Runge-Kutta 4 t,l -order method. 

In the laboratory coordinate system, the wave function is also defined by0 
in which it is necessary to replace amplitudes a m with v k m (z = h) from (101. 
Account should be taken of the fulfillment of the condition qm{q m — 2 ky) < Tg , 
which excludes exponentially attenuated (evanescent) waves and imposes certain 
restrictions on the admissible values of diffraction orders m and grating velocity 
V. In practice, diffraction of neutrons on a moving grating takes place when 
they pass through a rotating disk m that has radial rectangular grooves with 
a period d and depth h. If R is the distance of a ring with grating structure 
from the center of the disk and fd is its rotation frequency, then V = 2 irRfd- 

If the Bragg condition a\ = 0 is fulfilled and all orders except m = 0 and 
to = 1 can be ignored, then ( |10[ ) is reduced to the well-known X-ray crystal 
optics Takagi equations in the Laue geometry, whose solution has the following 
simple form: 'Fo(z) = A0cos(7rz/L ex ), 4 / i(z) = —iA0sin(Trz/L ex ). Here L ex = 
n2ko z /k 2 is the so-called extinction length at which the radiation from the zero- 
order fully pumped into the first order and vice versa. Under conditions of the 
experiments [2- 7 , it is a fraction of a micrometer which is comparable with the 
depth of the grating grooves h = 0.14 /tin. 

Figure [3] shows the intensities of the waves of three diffraction orders I m = 
|a m | 2 , where a m is defined in ([3]) (model 1), and J m = |'F m (/t)| 2 is defined within 
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the dynamical theory (model 2) in relation to the parameter c = 2hV/cLVo. It 
can be clearly seen that the intensity of the 0 th order in model 2 rises with 
increasing c (i.e. grating velocity) more rapidly than in model 1. The intensity 
of minus 1 st order in model 2 decreases faster than in model 1, which is in better 
agreement with experiments. We used the following data [5]: Vo = 4.52 m/s, h 
= 0.14 mm, d = 5 /im, R = 6 cm. At maximum rotation frequency of the disk 
fd =105 Hz the velocity V = 39.6 m/s and c max = 0.49. 



Figure 3: Intensities I rrl (solid curves, model 1) and J m (dashed curves, model 
2 ). 


Figure [4] shows the energy spectra 


1(E) 



\^ m (h)\ 2 exp[- 


(E - E m ) 2 
2<r 2 


( 11 ) 


calculated at different rotation frequencies fd, where the dispersion a 1 = 
a m + °a is determined by the convolution of the spectral functions of monochro¬ 
mator and analyzer of the gravitational spectrometer [2|, in which the scanning 
of energy E = mg AH is performed by vertical movement of the analyzer at a 
distance AH, where g is free fall acceleration. The spectral width of the vertical 
velocities AV/V 0 «0.02 |2j|4]. 

Above we assumed the depth of the grooves to be fixed (h = 0.14 /im). It is 
interesting to see how the intensities of the diffraction orders J m depend on h. 
As is shown on Fig. [5j for example, at a depth of h = 0.19 /im the zero-order 
intensity is fully suppressed, and the second-order intensity increases. At h = 
0.23 n m the intensities of minus 1 st and minus 2 nd orders are equal, and with 
further increasing h the minus 2 nd order becomes dominant. Let us note that 
J m (h) also depends on the frequency fd- 

Considering UCN diffraction on a grating, we assumed that a plane (i.e. 
non-divergent) wave falls on it. In fact, in typical experimental conditions [2j|5] 
there is a fairly broad (compared to the vertical velocity Vo) distribution of 
the neutron horizontal velocities v x . Its absolute value can be as high as 4.5 
m/s, which is determined by the value of the boundary velocity for the reflec¬ 
tions from the walls of the vertical neutron guide [3;. To take into account how 
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Energy (E - E 0 ), neV 

Figure 4: Energy spectra after diffraction at grating rotation frequencies fd= 40, 
75 and 105 Hz. 



this distribution affects on diffraction, it is necessary to replace in the above 
arguments ky with ky + q v , where q v = mv x /H - random increments of the 
x-component of the wave vector, and convolve d-shaped angular spectrum F(q) 
in © with the function of the velocity distribution. A specific feature of the 
diffraction on a moving grating is a decrease of angular divergence Ad of the 
incident beam in the grating coordinate system with an increase of the speed of 
its movement. This improves the condition A 9 < X/d for more reliable observa¬ 
tion of diffraction. The distribution of horizontal velocities v x with dispersion 
<r„ leads to an additional energy broadening a Em — k\q m \a v of diffraction or¬ 
ders. This broadening is independent on the grating velocity and grows with 
increasing order number. The typical values of (Jei ~l-2 neV. 
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